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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Si) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  con- 
verted to  metric  (Si)  units  as  follows: 

Multiply By To  Obtain 

degrees  (angle) 
inches 

kips  (force) 

pounds  (mass)  per  cubic  foot 


0.017 1+5329  radians 
2.5I+  centimetres 

1+U18.222  newtons 

16.018U6  kilograms  per  cubic  metre 


2 


CONTENTS 


Page 


PREFACE 1 

CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Si)  UNITS  OF 

MEASUREMENT 2 

CHAPTER  1 INTRODUCTION 5 

1.1  Background 5 

1.2  Objective 7 

1.3  Scope 7 

CHAPTER  2 PROBLEM  GEOMETRY,  BOUNDARY  CONDITIONS,  AND 

CONSTITUTIVE  AND  FIELD  EQUATIONS 8 

2.1  General 8 

2.2  Problem  Geometry 8 

2.3  Constitutive  Equations 9 

2.4  Transformation  of  the  Elastic  Properties 10 

2.5  Field  Equations 12 

2.6  Representation  of  the  Boundary  Condition l4 

2.7  Stresses  and  Displacements  in  Cylindrical 

Coordinate  Systems l4 

CHAPTER  3 DERIVATION  OF  THE  GENERAL  SOLUTION  OF  THE  PROBLEM 21 

3.1  General 21 

3-2  Stress  Function 22 

3.3  Complex  Representation  of  Stresses  and  Displacements 27 

3.4  Determination  of  the  Stress  Function 30 

3-5  Spatial  Stress  Distribution  for  Sample  Problem 34 

3.6  Spatial  Displacement  Distribution  for  Sample  Problem 35 

CHAPTER  4 CONCLUSIONS  AND  RECOMMENDATIONS 44 

REFERENCES 45 

APPENDIX  A:  NOTATION 47 


2.1  Problem  geometry 19 

2.2  Two-dimensional  element  from  the  boundary 20 

3.1  Distribution  of  radial  and  tangential  stresses  along  the 
radius  for  0=0  degrees  at  an  angle  of  inclination 

of  30  degrees 36 

3.2  Contour  for  radial  stress  at  an  angle  of  inclination 

of  30  degrees 37 

3.3  Distribution  of  tangential  stress  along  the  boundary  of 

the  borehole  whose  angle  of  inclination  is  30  degrees 38 

3.^  Distribution  of  radial  and  shear  stresses  along  the 
boundary  of  the  borehole  whose  angle  of  inclination 

is  30  degrees 39 

3-5  Distribution  of  radial  displacement  ur  along  the 
radius  for  0=0  and  90  degrees  at  an  angle  of 
inclination  of  30  degrees L0 

3.6  Contour  for  radial  displacement  at  an  angle  of 

inclination  of  30  degrees ll 

3.7  Distribution  of  radial  and  tangential  displacements 

along  the  boundary  of  the  borehole  whose  angle  of 
inclination  is  30  degrees U2 

3.8  Contour  for  tangential  displacement  at  an  angle  of 

inclination  of  30  degrees b3 


LIST  OF  TABLES 

Table  Page 

2.1  Direction  cosines l6 

2.2  Values  of  q-y  in  the  formulas  of  transformations 

(Equation  2.3) l6 

2.3  Value  of  C^j  from  Equation  2.1 17 

2.b  Value  of  C-y  from  Equation  2.3 18 


RESPONSE  OF  LINEAR  ELASTIC  TRANSVERSE- ISOTROPIC  MEDIA 
TO  BOREHOLE  PRESSUEEMETER  LOADINGS 

CHAPTER  1 
INTRODUCTION 

1.1  BACKGROUND 

To  predict  the  ground  shock  from  surface  or  aboveground  nuclear 
detonation,  research  on  the  use  of  two-dimensional  finite  difference 
wave  propagation  codes  that  treat  nonlinear  hysteretic  media  is  being 
conducted.  To  use  these  codes,  the  constitutive  properties  (stress- 
strain  and  strength)  of  the  in  situ  earth  media  must  be  determined  for 
fast  loading  rates  in  the  unconsolidated-undrained  state.  Conven- 
tionally, this  is  done  by  obtaining  undisturbed  samples  from  the  site 
and  testing  them  in  the  laboratory.  Inevitably,  the  in  situ  properties 
are  altered  to  some  extent  by  the  sampling  process.  Hence,  using 
in  situ  field  tests  that  give  some  indication  of  the  in  situ  constitutive 
properties  is  desirable. 

The  borehole  pressuremeter  (References  1-3),  which  measures  the 
increase  in  volume  per  unit  length, AV  , of  a borehole  under  an  in- 
creasing uniform  internal  pressure,  Pq  , is  one  of  several  tests  that 

can  be  used  to  infer  information  about  the  in  situ  constitutive  orop- 
2 

erties.  This  device  has  been  used  in  the  constitutive  property  in- 
vestigation for  several  high-explosives  (HE)  tests  (References  4-7)  used 
as  test  cases  to  study  the  accuracy  of  and  the  necessary  improvements  in 
the  ground  shock  prediction  procedure. 

In  homogeneous  isotropic  linear  elastic  materials,  the  shear  modu- 
lus G can  be  directly  determined  from  the  borehole  pressuremeter  test: 


Symbols  used  in  this  report  are  listed  and  defined  in  the  Notation 
2 (Appendix  A). 

As  with  any  test,  it  has  its  own  sources  of  error,  the  most  prominent 
of  these  being  volume  measurement  errors  leading  to  alteration  of 
the  constitutive  properties  due  to  the  drilling  of  the  borehole. 
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V P 

G = ir  (1*1} 

where  V is  the  initial  volume  of  a unit  length  of  the  borehole, 
o 

However,  real  earth  materials  are  often  highly  anisotropic.  The  inter- 
pretation of  data  from  J.aboratory  and/or  field  teste  based  on  a mathe- 
matical constitutive  relation  that  does  not  account  for  anisotropy  may 
lead  to  erroneous  conclusions. 

The  most  common  departure  from  a state  of  isotropy  in  an  earth  ma- 
terial is  layering  or  stratification  during  its  deposition.  This  is  the 
case  whether  induced  by  natural  causes,  such  as  sedimentary  deposits,  or 
in  the  construction  of  fills  where  the  earth  materials  are  placed  and 
compacted  in  horizontal  lifts.  For  these  conditions,  although  marked 
differences  may  be  noted  between  the  vertical  and  the  horizontal  direc- 
tions, generally  no  direction  preference  will  exist  in  the  horizontal 
planes.  Such  a material  is  said  to  be  transversely  isotropic 
(References  8 and  9). 

For  the  linear  elastic  transverse-isotropic  material,  five  material 
constants  are  needed  to  completely  describe  material  behavior  (Refer- 
ences 8 and  9).  A technique  for  the  determination  of  some  or  all  of 
these  five  material  constants  in  the  field  is  of  great  importance  to  the 
material  properties  investigator.  Such  a technique  (Reference  10)  would 
give  an  early  deduction  of  those  sites  that  are  strongly  anisotropic  and 
provide  data  for  use  in  the  fitting  of  a transverse-isotropic  model 
(Reference  7)  for  the  materials.  It  appears  that  a series  of  pressure- 
meter  tests  in  boreholes  inclined  at  several  different  angles  to  the 
axis  of  symmetry  of  the  material  will  yield  this  information.  Hence, 
there  was  a need  to  obtain  an  analytical  solution  for  the  inclined  bore- 
hole pressuremeter  problem.  Only  the  special  case  of  this  problem  for 
a transverse-isotropic  material  in  which  the  axis  of  the  test  (i.e.,  the 
axis  of  a pressurized  cylindrical  cavity  of  infinite  extent)  is 
perpendicular  to  the  plane  of  isotropy  has  been  solved  analytically 
(Reference  9). 
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1.2  OBJECTIVE 


The  objective  of  this  investigation  is  to  develop  a general  closed 
form  solution  to  the  borehole  pressuremeter  problem  where  the  axis  of 
the  borehole  is  inclined  to  the  plane  of  isotropy  of  an  elastic 
transverse-isotropic  medium.  Such  a solution  can  he  used  to  determine 
stresses  and  displacements  in  the  medium  in  terms  of  the  transverse- 
isotropic  properties  of  the  material.  It  is  possible,  therefore,  to  de- 
duce the  transverse-isotropic  properties  of  the  materials  in  terms  of 
the  volume  change  of  the  borehole. 

1.3  SCOPE 

The  problem  geometry , boundary  conditions , and 
field  equations  are  presented  in  Chapter  2.  Chapter 
ivation  of  the  general  solution  of  the  problem  and  a 
Conclusions  and  recommendations  are  given  in  Chapter 


constitutive  and 
3 contains  the  der- 
s ample  problem. 

k. 


CHAPTER  2 


1 


PROBLEM  GEOMETRY,  BOUNDARY  CONDITIONS,  AND 
CONSTITUTIVE  AND  FIELD  EQUATIONS 

2.1  GENERAL 

Because  the  properties  of  the  material  are  directionally  dependent , 
it  is  more  conven'c-nt  to  obtain  the  solution  of  the  problem  in  a 
Cartesian  coordinate  system  first  and  then  through  a coordinate  trans- 
formation to  determine  the  stresses  and  the  displacements  in  a cylin- 
drical coordinate  system.  Therefore,  three  coordinate  systems  are 
needed.  The  first  is  a Cartesian  coordinate  x'y'z'  in  which  x'y' 
is  parallel  to  the  plane  of  isotropy  (Figure  2.1)  and  for  which  the 
constitutive  equations  for  a linear  elastic  transverse-isotropic  material 
are  well  known  (References  8 and  9).  The  second  coordinate  system  is 
also  Cartesian,  xyz  , in  which  the  solution  of  the  borehole  pressure- 
meter  problem  is  obtained  (Figure  2.1).  Finally,  because  the  problem 
is  axisymmetric , it  is  convenient  to  transform  the  final  results  from 
the  xyz  coordinate  system  to  a cylindrical  coordinate  system,  r0z 
(Figure  2.1b). 

2.2  PROBLEM  GEOMETRY 

The  geometry  of  the  problem  is  shown  schematically  in  Figure  2.1. 
Figure  2.1a  shows  a three-dimensional  view  of  the  problem  with  the  rela- 
tive position  of  the  Cartesian  coordinate  systems  x’y'z'  and  xyz  . 
Figure  2.1b  shows  a two-dimensional  detailed  view  of  the  problem  relative 
to  both  the  cylindrical  coordinates  rOz  and  the  Cartesian  coordinates 
xyz  . The  axis  of  symmetry  of  the  material  and  the  axis  of  symmetry  of 
the  cylindrical  cavity  are  assumed  to  intersect  at  an  angle,  0 (Fig- 
ure 2.1a).  Therefore,  the  intersection  of  the  cylindrical  cavity  with 
the  plane  of  isotropy  forms  an  ellipse  A'BC'D  in  the  x'y'  plane  (Fig- 
ures 2.1a  and  2.1b).  The  intersection  of  the  cylindrical  cavity  with  the 
r0-  or  xy-plane  is,  of  course,  a circle  ABCD  (Figure  2.1a).  Further- 
more, it  is  assumed  that  the  radius  of  the  cylindrical  cavity  is  b and 
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that  its  surface  is  under  normal  stress,  Pq  , which  does  not  vary  along 
the  cavity. 

2.3  CONSTITUTIVE  EQUATIONS 

Let  the  x'y’  plane  of  an  x'y'z'  coordinate  system  (Figure  2.1a) 
be  the  plane  of  isotropy  of  the  material.  The  constitutive  equations  for 
a linear  elastic  transverse-isotropic  material  (References  8 and  9)  are: 

_ _1  v v_|_  ^ 

£x'  “ E ax'  " E V ~ E'  V 


_ . 1 vj_ 

£y'  - ~ E °x'  E °y'  - E'  °z' 


(a  ,+  a S + — a , 
V x’  y ) E'  z’ 


(2.1) 


where 


y'z'  2G' 


'x'y'  " 2G 


a , , a , , a , = total  normal  stress  components  parallel  to 


x’  y'  z 


x'-,  y'-,  and  z'-axes,  respectively 


°x'z'’  ° 'z' 5 °x'  ' = shearing  stress  components  in  x'z'-, 

y'z'-,  and  x'y '-planes,  respectively 

, s ,,  e^,  = total  normal  strain  components  parallel  to 
x'-,  y'-,  and  z'-axes,  respectively 

ex'z’’  ev'z'’  Gx'v'  = ’bo"tal  shearing  strain  components  in  x'z'-, 

y'z'-,  and  x'y '-planes,  respectively 

E = Young's  modulus  in  the  plane  of  isotropy 

E'  = Young's  modulus  in  a plane  normal  to  the 
plane  of  isotropy 
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v = Poisson's  ratio  that  characterizes  the  trans- 
verse reduction  in  the  plane  of  isotropy  due 
to  stress  in  the  same  plane 

v'  = Poisson's  ratio  that  characterizes  the  trans- 
verse reduction  in  the  plane  of  isotropy  due 
to  stress  normal  to  it 

G'  = shear  modulus  for  a plane  normal  to  the  plane 
of  isotropy 

g 

G = ~rn — ; r = shear  modulus  for  the  tilane  of 

2(1  + v) 

isotropy 

The  elastic  properties  that  appear  in  Equations  2.1  depend  on  the 
direction  of  the  axes  of  the  chosen  coordinate  system.  If  the  direction 
of  the  axes  varies,  then  the  elastic  properties  vary.  Only  in  the  case 
of  an  isotropic  body  the  elastic  properties  are  invariant  in  any  orthog- 
onal coordinate  system.  However,  there  are  always  unique  relationships 
of  the  elastic  properties  in  one  coordinate  system  to  the  elastic  prop- 
erties in  another  coordinate  system.  These  relarionships  could  be  de- 
rived through  transformation  formulas  that  transform  one  coordinate 
system  into  another.  Therefore,  the  elastic  properties  that  appear  in 
Equations  2.1  for  the  coordinate  system  x'y'z'  could  be  transformed 
into  the  elastic  properties  for  the  coordinate  system  xyz  (Figure  2.1) 
through  transformation  formulas. 


2. U TRANSFORMATION  OF 

THE  ELASTIC  PROPERTIES 

Let  C. , be  the  elastic  properties  for  the  coordinate  system  xyz 
• 0 

and  let  0!^  be  the  elastic  properties  for  the  coordinate  system  x'y'z' 
(Figure  2.1a).  The  position  of  the  coordinate  system  xyz  with  respect 
to  the  coordinate  system  x'y'z'  is  defined  by  Table  2.1  and  the  follow- 
ing relations: 

x = x' 

y = y'  cos  + z'  sin  ^ (2.2) 

z = -y'  sin  rp  + z'  cos  ip 
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The  transformation  formulas  that  relate  C.  to  C' , 

x ij  iJ 

Reference  9 and  can  he  written  as 


are  given  in 


C = C'  q .q  , 
lj  mnTni  nj 


(2.3) 


The  values  of  q.  are  defined  in  Table  2.2  where  the  first  index,  i , 

indicates  the  number  of  the  row  and  the  second  index,  j , shows  the 

number  of  the  column.  Thus,  q denotes  the  element  belonging  to  the 

ith  row  and  column;  for  example,  q^  = 62  , q^  = ’ 

q,.^  = 0,<5O  + 0o<5-,  > and  so  forth.  The  values  of  6 , n and  0 

59  1221  n n n 

(n  = 1,  2,  3)  are  given  in  Table  2.1;  the  values  of  C!  can  be  obtained 

— J 


from  Equations  2.1  and  are  given  in  Table  2.3;  and  the  values  of  C. 


ij 

obtained  from  Equation  2.3  are  given  in  Table  2.1*. 

Having  determined  the  value  of  the  elastic  properties  C. . , the 

* J 

general  constitutive  equation  for  a linear  elastic  transverse-isotropic 
material  in  an  xyz  coordinate  system  may  be  written  as: 


1 /v'  2 v 2 \ 

:x  = E Cx  “ VF  Sin  * + I COS  V °y 

/v'  2 v 2 \ /v'  v \ 

- ^ cos  * + £ sin  ^ + (jr  ~ e)  sin  2* 

v * - (£ ai] 

/ 1 2v’\  2 2 ~\  f v'  , It  k . 

+ jET ) cos  4 sin  4j  oy  + ^7  (sin  4 + cos  4) 


2 v 2 
sin  4 + — cos  4 


yz 

1* 


’)  °x  + [e  Cos1+  * + E7  Sin"  ^ 


y (2.1*) 


a + 1__  u 

\E  E’  GV 


2 . 2 
cos  4/  sin  ip 


] * [(!  * w) 


2 , 
cos  4' 


/I  v'  \ 2 1 1 

Ve5-  + WJ  Sin  * " 2G^  C0S  2*J  Sin  2 * ayz 


NOTE:  Equations  2.1*  are  continued  on  following  page 


Indices  assume  values  1,2... 6.  A repeated  index  is  to  be  summed  over 
its  range.  Quantities  are  referred  to  rectangular  Cartesian  coordi- 
nates x.  . 
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sin  ip  + cos  ip 

k . . 1 U 


= - (jT  cos2  IP  + | sin^  ax  + [-  ~ ( 

+ (I  + F " f)  COs2  * Sin2  *]  °y  + [I  si 
+ (f  " F")  cos2  * sin2  *]  °z  + [(e  + f)  sin' 


k . . . u ,) 


in  ip  + ^7-  cos  ip 


*P 


2 1 1 
cos  ip  + 2qT  cos  20 J san  2\t>  a. 


’ I (f  ' I)  Sln  2*ax*2  [(I  * f) 

- (f-  + f)  «i 


2 , 
COS  Ip 


in2  ip  - cos  20  J sin  20 


I [(I  + f)  3in2  * " (f  + f) 


2 , 
COS  0 


+ 2F  COS  2*J  Sin  20  aZ  + 2G^  ayz 


1 2 1 • 2 . 
exz  = IF  COS  ^ + 2G  Sln  * 


1 . 2 . , 1 2 . 
exy  = 2F  Sin  * + 2G  C°S  * 


^ 2.  it 

(cont'd) 


where 


a 

z 


total  normal  stress  components  parallel  to 
x-,  y-,  and  z-axes,  respectively 


total  shearing  stress  components  in  xz- , yz-, 
and  xy-planes,  respectively 

total  normal  strain  components  parallel  to 
x-,  y-,  and  z-axes,  respectively 


E 

XZ 


yz 


E 

xy 


total  shearing  strain  components  in  xz- , yz-,  and 
xy-plenes,  respectively 


2.5  FIELD  EQUATIONS 

In  the  case  of  small  displacements  of  a continuous  body,  the  rela- 
tionships between  the  components  of  strain  and  displacements  (Refer- 
ence 9)  are: 
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r 


_ 3u  3v  3w 

x 3x  ’ ey  3y  ’ Ez  3z 

P - i ( iiL  + ,.  _ 1 /3v  3w  \ , . 

xz  2 \3x  3z/  ’ eyz  2 \3z  3y / > (2>5) 

= 1 /3u  3v\ 

xy  2 \3y  3x/  J 

where  u , v , and  w are  the  displacements  in  the  x-,  y-,  and  z- 

directions,  respectively. 

In  the  problem  under  consideration  (Figure  2.1),  the  stresses  and 
displacements  are  independent  of  z and  become  functions  of  x and  y 
alone.  Therefore,  Equations  2.5  can  be  written  as: 


g 3u(xjX).  =i vObil  . _n 

x 3x  ’ Gy  3y  ’ Ez  “ ° 


, _ 1 3w(x,y 

'xz  2 3x 


_ 1 3w(x,y ) 
"y z 2 3x 


= I 3v(x,y)~| 

2 L 3y  3x  J 


(2.6) 


Equations  2.6  leads  to  the  compatibility  equations  that  guarantee  the 
body  is  continuous. 

The  stress  components  in  a continuous  body  in  equilibrium  Tinder  the 
action  of  surface  and  body  forces  satisfy  three  differential  equations 
of  equilibrium.  In  the  case  under  consideration,  these  equations  take 
the  following  form: 


3 a do 

+ —M.  = o 

3x  3y 


3ct  3ct 

—352.  + _2L  = o 

3x  3y 


(2.7) 


2.6  REPRESENTATION  OF 

THE  BOUNDARY  CONDITION 

Let  X and  Y Be  the  x and  y components,  respectively,  of  a 
distributed  surface  force  per  unit  area;  the  boundary  stress  equations 
(Figure  2.2)  can  be  written  as  (Reference  9): 


X = 

0 

X 

cos 

(n,x) 

+ 0 

xy 

COS 

(n,y)  + °xz 

cos 

(n,z) 

Y = 

Q 

X 

*< 

cos 

(n,x) 

+ 0 

y 

cos 

(n,y)  + oyz 

cos 

( n , z ) 

0 = 

0 

xz 

cos 

(n,x) 

+ 0 

yz 

cos 

(n,y)  + oz 

cos 

(n,z) 

For  the  above  equations,  the  following  relationships  exist: 

cos  ( n , x ) = 

ds 

i \ dx  V 

cos  (n,y)  = - — 

cos  (n,z)  = 0 


(2.8) 


(2.9) 


2.7  STRESSES  AND  DISPLACEMENTS  IN 
CYLINDRICAL  COORDINATE  SYSTEMS 

The  relations  between  the  stresses  and  the  displacements  in  the 
Cartesian  and  cylindrical  coordiante  systems  with  the  same  z-axis  (Fig- 
ure 2.1)  are: 


2 2 

o =o  cos  0+o  sin  0 + 2o  cos  0 sin  0 

r x y xy 

2 2 

o_  = o sin  0+o  cos  0 - 2o  cos  0 sin  0 

0 x y xy 

o =o 
z z 

o = (a  - o \ cos  0 sin  0 + o (cos2  0 - sin2  0) 
r0  \ y x)  xy  v ' 


> (2.10) 


o =o  cos  0 + o sin  0 
rz  xz  yz 


o„  = - a sin  0 + a cos  0 
0z  xz  yz 


1 


and 


where 


u = u cos  0 + v sin  0 
r 


vQ  = - u sin  0 + v cos  0 


(2.11) 


0 = tan 


-1  Z 


(2.12) 


and  ur  and  vQ  are  the  radial  and  tangential  displacements,  respec- 
tively. Therefore,  if  the  stresses  and  the  displacements  in  the  xyz 
coordinate  system  are  known,  the  corresponding  stresses  and  displace- 
ments in  the  r0z  coordinate  system  can  easily  he  obtained. 
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A 


Table  2.2.  Values  of  q^j  in  the  formulas  of 
transformations  (Equation  2.3). 


X 

1 

2 

3 

h 

5 

6 

1 

61 

4 

4 

26263 

2636i 

26 16  2 

2 

nl 

4 

n3 

2T]2T)3 

2n3nl 

2nln2 

3 

4 

02 

4 

20203 

20391 

201G2 

U 

nl0l 

n202 

n303 

n203  + 

n302 

nl03  + n30l 

r'l02  + 

n2Gl 

5 

Vi 

®262 

V3 

V3  + 

0362 

Gl63  + G36l 

Gl62  + 

Vl 

6 

6lnl 

62n2 

53n3 

62n3  + 

63T,2 

6ln3  + 63nl 

6ln2  + 

62n! 

L 


OlM 


Table  2.4.  Value  of  C. , from  Equation  2.3. 


CHAPTER  3 

DERIVATION  OF  THE  GENERAL  SOLUTION  OF  THE  PROBLEM 

3-1  GENERAL 

The  problem  of  the  determination  of  stresses  and  displacements 
around  an  infinite  cylindrical  cavity  can  he  formulated  analytically  by 
use  of  the  equations  presented  in  Chapter  2.  The  solution  will  be 
unique  if  it  satisfies  the  equilibrium  and  compatibility  equations  and 
the  boundary  conditions  (Reference  11).  This  is  done  below  by  expressing 
stresses  and  displacements  as  complex  harmonic  functions  in  a manner 
similar  to  that  developed  by  Lekhnitskii  (Reference  9).  Recalling 
Equations  2.4, 


e = C,  ..a  + C,o0  + C,  _cr  + C-.,cr 

x 11  x 12  y 13  z 14  yz 


e = c,  a + C a + C„_c  + C~i  o 

y 12  x 22  y 23  z 24  yz 


e = C,  _a  + C„_a  + C o + C_,a 
z 13  X 23  y 33  Z 34  yz 


(3.1) 


2e  = C,,  a + Co)  a + Co)  a + Cua 
yz  14  x 24  y 34  z 44  yz 


2e  = Ccca 
xz  55  xz 


2e  = C^a 
xy  66  xy 


where  the  values  of  C are  given  in  Table  2.3.  Since  e =0  (see 

i (3  ^ 

Equations  2.6),  the  third  equation  of  3.1  leads  to 


a = - 7—  (C,~o  + Coo0  + C,,  a ^ 
z C33  V 13  x 23  y 34  yz^ 


(3.2) 


Substitution  of  Equation  3-2  into  Equations  3.1,  gives 
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in  which 


e = a,  ,0  + a,  „o  + a,  < o 

x 11  x 12  y 14  yz 


e = a,  _o  + aoo0  + a„,  o 
y 12  x 22  y 24  yz 


2e  = a,  .0  + aol  o + aM  o 

yz  14  x 24  y 44  yz 


2e  = acco 
xz  55  xz 


2e  = a,-  ro 
xy  66  xy 


a - r a3  jL3. 
ij  iJ  C33 


Equations  3.3  can  he  written  in  terms  of  the  displacements  as 

3u 

— = a,  .,0  + aio0 

3x  11  x 12  y 

3v 

r — = a.  o0  + aoo0 
3y  12  x 22  y 

3w 

t — = a,  i o + a_i  ,0 
3y  14  x 24  y 

3w 

. _0 
3x  55  xz 

3u  3v  _ 

3y  + 3x  a66°xy 

3.2  STRESS  FUNCTION 

Equations  2.7  (the  equilibrium  equations)  can  be  satisfied  for  a 
homogeneous  medium  by  the  introduction  of  the  following  stress  functions 
(References  9»  12,  and  13): 


+ a, , o > 
14  yz 


+ aoi  0 

24  yz 


+ a.  , o r 
44  yz 


(3.5) 
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The  compatibility  equations  can  be  satisfied  by  substitution  of 


Equations  3.6  into  Equations  3.5  and  elimination  of  u , v , and  w by 
differentiation  (Reference  13).  Therefore,  the  following  system  of  dif- 
ferential equations  that  the  stress  functions  must  satisfy  can  be  easily 
obtained: 


L,P(x,y)  + L,Q(x,y)  = 0 


L^PU.y)  + L2Q(x,y)  = 0 


(3.7) 


where  , an<^  ^ are  dif’f’erential  operators  of  the  second, 

third,  and  fourth  orders,  respectively,  that  have  the  form: 


2 2 
A 

L2  = aUU^2+a55^T 


93  93 

L - - a.  ^ - a ^ 

3 9xJ  4 9x9y^ 


Lb  = a22  + (2ai2  + a66)  zty  + ail  ZJ 


_9_ 

3x29y‘ 


9_ 
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(3.8) 


J 


For  the  components  of  stresses  and  displacements  around  the  cylin- 
drical cavity  to  be  continuous  and  single-valued  functions  of  the 
coordinates  xyz  , the  stress  functions  P(x,y)  and  Q(x,y)  must  sat- 
isfy Equations  3.7  and  the  boundary  conditions. 

The  general  differential  equations  in  terms  of  P(x,y)  and 
Q(x,y)  , separately,  can  be  obtained  by  application  of  the  operator  L? 
on  the  first  equation  of  the  system  3.7  and  the  operator  L_  on  the  sec 
ond  equation  and  subtraction  of  the  results.  Thus: 


(LiL2  ~ Ll)  p(x,y)  = 0 (3.9) 

Similarly, 

(lul2  ~ ^3}  Q(x>y)  = 0 (3.10) 


Equations  3.9  and  3.10  are  sixth  order  differential  equations  where  the 
operator  of  the  sixth  order  L^L2  ~ can  ^-ecomPose<i  into  six 
linear  operators  of  the  first  order.  Hence,  Equations  3.9  and  3.10  can 
be  represented  in  the  following  forms: 


and 


in  which 


D6D5D4D3D2D1  P(x’y)  = 0 
D6D5DUD3D2D1  Q(x’y)  = 0 


J 


(3.11) 


3_ 

3y 


3 

yk  3x 


(k 


1,2. ..6) 


(3.12) 


where  represents  the  roots  of  the  following  algebraic  equation  that 

corresponds  to  the  differential  Equations  3.9  and  3.10: 


Lj^p)  L2(p)  - L^( P ) = 0 


(3.13) 


According  to  Equations  3.8,  L^p)  , L^(p)  , and  L^(p)  can  be  written 

as : 


2h 


(3.114) 


L2(y)  = a^p2  + %k 

L3(y)  = -alUy2  - a^  \ 

Lu(y)  = aup^  + (2a12  + ag6)  y2  + a22  ^ 

Three  of  the  roots  of  Equation  3.13  are  independent;  the  other  three  are 
their  complex  conjugates  (Reference  9). 

The  integration  of  Equations  3-9  and  3.10,  therefore,  can  be  reduced 
through  Equations  3.11  to  the  integration  of  six  equations  of  the  first 
order.  The  general  integral  is  equal  to  functions  of  the  arguments 

\ = x + vn 

and  \ (k  = 1,  2,  3)  (3.15) 

\ = x + V J 


and  can  be  written  as 


QU,y)  = Ji  [M^)  + ] 


(3.16) 


where  y^  is  the  complex  conjugate  of  y^  , 
jugate  of  2^  , and  P^z^J  and  Q^z^J  are 

Pk(Zk)  and  Sc(Zk)  ’ resPectively- 

Since  the  functions  P(x,y)  and  Q(x,y) 
3.8,  the  following  relations  between  P(x,y) 


z,  is  the  Comdex  con- 
k 

the  complex  conjugates  of 

satisfy  Equations  3.7  and 
and  Q(x,y)  exist: 


MM ' - 


ap>(\) , 

dz. 


Vi 


+ b. 


(3.17) 
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or 


"kW 


+ Bk 


(3.18) 


where  a^  , b^  , , and  are  arbitrary  constants.  Hence,  the 

stress  functions  P(x,y)  and  Q(x,y)  (Equations  3-l6)  can  be  written 
as : 


P(x,y)  - P^)  . * P2(z2)  * F^} 

+ P3(b)+  W 


(3.19) 


Q(x,y)  = ) 


X1P1' 


l(Zl)  + X2P2  (Zp) 


X2P2 


where 


and 


V 


+ b. 


(3.20) 
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Therefore,  the  general  solution  to  the  borehole  pressuremeter  problem 

can  be  completely  determined  by  determining  the  functions  P (z  ) . But 

k k 

before  this  can  be  done,  the  stresses  and  the  displacements  have  to  be 

expressed  as  functions  of  P (z,  ) . 

K.  K 

3.3  COMPLEX  REPRESENTATION  OF 
STRESSES  AND  DISPLACEMENTS 

Since  the  stresses  are  functions  of  the  second  derivative  of  P(x,y) 
(Equations  3.6),  and  the  displacements  are  functions  of  the  first  deriv- 
ative of  P(x,y)  (Equations  3-5)»  it  is  more  convenient  to  introduce  the 
new  functions  of  the  complex  variable  z 


dP, 


£*)  = p. 


(\)  = -ir^  = pk(zk) (k  = i* 2) 


1 dP3(Z3)  1 

|(Z3)=^^  =i;p3(z3) 


(3.22) 


With  the  help  of  these  functions,  the  expressions  for  the  first  and  sec- 
ond derivatives  of  P(x,y)  and  for  the  firBt  derivatives  of  Q(x,y) 
with  respect  to  x and  y may  be  written  in  the  following  way: 


ff.(x,y)  = 

3x 


MZl)  + *l(Zl)  + ^2  (Z2 ) + MZ2) 


+ V3(Z3)  + SWs)  (3-23) 


ll^LL 

8x2 


* *i(2i)  * *i(z2)  * Wd 


l3*iN  + 


3 3 


3v3/ 


3P. 


k(zk)  dPk(zk)  dzk  ,,  \ / . 

3y  dz  dy  ~ Pk  k(Zk)  ~ Mk\(Zk) 
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(3.21+) 


(3.25) 


3P(x ,y ) = 

3y 


32p(*>y) 

..2 


ui*i(zi) 

"1*1 (zi) 


* '•'l*l(zl)  + "2*2 (z2)  + "2*2(z2) 

* >,3»3+3(z3)  * 

+ IV>;  (zi)  * 4*2 (Zs)  * 4*s(zs) 

* u3l3»3(z3)  * 


(3.26) 


(3.27) 


3 P(x,y) 


3x3y  = yl*i(Zl)  + Ml*i(zl)  + + y2^(Z2) 


+ y3X3^(Z3)  + y3V3(Z3)  (3'28) 


Q(x,y)  = X1*1(zJ  + + A^2(z2)  + 


+ ^3  (Z3)  + S(Z3)  (3'29) 


3Q(x,y) 

9x 


- vi(zi)  * vIFIJ + * vJF^y 


*Kzs)  * (3-30) 


3Q(x,y) 

3y 


ylAl^(Zl)  + Wl)  + x2v2^(*2)+ 


X2V2*2(Z2) 


+ V3(Z3)  + Vlh)  (3-3l) 


According  to  Equations  3.2  and  3.6,  and  on  the  "basis  of  Equa- 
tions 3.23  through  3.31,  the  general  expression  for  the  components  of 
stresses  can  he  obtained  as: 
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°x  = 2Real  [ul*i(zl)  + 4^2)  + 

°y  = 2Real  [*i(zl)  + ^2 (z2 ) + X34,3(Z3)] 

°xy  = -2Real  [Vl(zl)  + v2+2(z2)  + V3*3(Z3)] 

°xz  = 2Real  |>lXl*i(Zl)  + y2X2^(Z2)  + Vsfc)] 
ayz  = -2Real  [A-^zJ  + ^(*2)  + *^3)] 


a = - — — (c.-0  + C00o  + CL,  a ^ 

z C \ 13  x 23  y 34  yz/ 


(3.32) 

(3.33) 
(3.3U) 

(3.35) 

(3.36) 

(3.37) 


The  displacements  u , v , and  w can  be  obtained  by  the  substitu- 
tion of  Equations  3-32  through  3.37  into  Equations  3.5  and  integration 
of  the  resulting  equations.  Thus: 

k=2 

u = 2Real  k?i  (*llMk + *12 " Xk&1^  *k(Zk) 

+ 2Real{[X3(allU3  + al2)  - alJ  ♦sft)}  (3-38) 


k=2  / a A \ 

v = 2Real  £ ^ + — - - a^j  *k(zk) 


+ 2Real  < X3  a12y3  + T 


a2l*  Xk 


2k  w \ 
TJ  *3(z3) 


(3.39) 


v =2Real  £ ^ *hh)  ♦kfk) 

(r  / a2i*\  auu' 

+ 2Real  + “Sj)  - — j ♦sfs) 


(3. Ho) 


The  stresses  and  the  displacements  in  cylindrical  coordinate  systems 
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can  be  obtained  by  substitution  of  Equations  3.32  through  3.40  into 
Equations  2.10  and  2.11. 

It  is  clear  from  Equations  3.32  through  3.40  that  ) is  the 

only  function  needed  for  the  determination  of  stresses  and  displacements. 
This  function  can  be  determined  from  the  boundary  conditions,  as  shown 
in  the  next  section. 

3.U  DETERMINATION  OF 

THE  STRESS  FUNCTION 

The  relationships  between  the  stresses  along  the  boundary  and 
inside  the  region  can  be  obtained  by  the  combination  of  Equations  2.8 
and  Equations  3.32  through  3.37*  Thus: 

Mzi) + *i(zi) + Mz2)  + *2(z2) + VsN 

s 

+ A3$3(z3)  = J - Yds  (3.41) 

0 

yl*l  (Zl)  + Wl^l(Zl)  + v2*2(z2)  + w2*2(z2)  + V3*3(Z3) 

S 

+ p3X34»3(z3)  = J XdS  (3.42) 

0 


Al*l(2l)  * Vl(2l)  * *2*2  (z2)  * V2U2)  * *3 (Z3> 

**3H)'Cc  (3'1,3) 

where  S is  an  arc  length  along  the  boundary  and  Cq  is  an  arbitrary 
constant . 

The  arguments,  z,  , in  the  above  functions  can  be  written  as  (see 
K 

Equation  3.15): 

\ = j (1  - iuk)  exp(i0)  + "I  (1  + iwk)  exp(-iO)  (k  = 1,2,3)  (3.44) 
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A comparison  between  Equations  3.41  through  3.43  and  3.45  through  3.47 
leads  to 


S 0 « / . „ 

r ~ r _ v / knB  — ;in 

J J ^ Vn  bn  n bn 

0 0 n=l 


in0  -in0 

" — + b — 

bn  n br 


C C)  oo 

r _ r _ K?  / ^n®  -in0' 

/ XdS  = b / Xd0  = > [d  - — + d 5 — 

J J z.  In  n n n 

0 0 n=l  ' 


(3.48) 


(3.49) 


The  coefficients  bn  and  d^  car.  be  obtained  by  use  of  the  properties 
of  the  Fourier  series  (Reference  13)  and  the  problem  geometry  defined  in 
Section  2.2,  which  yield  the  following  values: 


J 


Substitution  of  Equations  3.50  into  Equations  3.^+5  through  3.^7  leads  to: 


Reai  px(2l)  + *2(z2)  + V3(z3)]  = Real  \"  “§7~  *7  (3'51) 

X"1  + “2*2(z2)  + “3i3*3(z32  = E'al  (3'52> 


Beal  [v.jfZj)  * + «,(*,)]  = 0 (3.53) 

When  Equations  3*51  through  3-53  are  solved,  this  results: 


P b 

■iW*  i 


2 r 


2r 


("l  - va) : (ll  - 1 

U1  * V3(“l  - “3)  * - *2) 


(sin  Q + j cos  Q 
sin  O - Ug  cos  0 


(3.58) 


3(Z3)  = 


P b 
o 


2 r 


2r 


(wg>l  - u >a)  +(»  - >,)  i 

"1  * V3(U1  - "3)  * VsC-s 


)J 


/ sin 

0 + 

i cos 

0 ' 

l sin 

0 - 

p^cos 

0 

> 

(3.59) 


The  distribution  of  stresses  can  be  determined  from  Equations  2.10,  3.32 
through  3.37  and  3.57  through  3.59»  and  the  distribution  of  displacements 
can  be  determined  from  Equations  2.11,  3.38  through  3.*+0,  and  3.5*+ 
through  3.56.  The  computer  program  BOREHOLE  was  developed  to  solve 
numerically  the  above  system  of  equations  and  to  generate  various  plots 
of  stress  and  displacement  distributions  around  the  cylindrical  cavity. 
Examples  of  the  distribution  of  stresses  and  displacements  are  given  in 
Figures  3.1  through  3.8. 

The  volume  change  of  a unit  length  along  the  generator  of  the  bore- 
hole can  be  obtained  from  the  radial  displacement  (Equation  3-38)  at 
r = b , and  0=0,  and  0 = n/2  : 

Ay  ir[b  + u(b,0)]  [b  + u(b,  j)]  - irb^ 

~ = ,2 

irb 

or 


AV 

V 


u(b,0)  + u(b,  u(b,0)  u(b, 

b + b2 


(3.60) 


Equation  3.60  is  a function  of  the  five  material 
well  as  the  angle  of  inclination  of  the  borehole,  ij; 


properties  as 
(Figure  2.1) . 
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Therefore,  the  solution  of  Equation  3.60  for  the  material  properties  is 
not  straightforward  and  requires  an  iterative  scheme  and  a large  computer 
program  such  as  BOREHOLE.  However,  the  solution  is  relatively  simple  if 
four  material  properties  as  well  as  the  volume  change  are  known. 

In  the  following  section,  spatial  stress  and  displacement  distribu- 
tions for  a sample  problem  are  investigated.  The  material  properties 
used  in  this  sample  problem  as  well  as  the  angle  of  inclination  of  the 
borehole  are  tabulated  below. 
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3 . 5 SPATIAL  STRESS 

DISTRIBUTION  FOR  SAMPLE  PROBLEM 

Figure  3.1  shows  a typical  result  of  the  radial  and  tangential 

stresses  along  the  radius  for  0=0  degrees'1’  (Figure  2.l)  at  an  angle 

of  inclination  of  30  degrees.  The  solid  line  shows  the  radial  stress 

while  the  dashed  line  shows  the  tangential  stress. 

Figure  3.2  shows  a radial  stress  contour  in  dimensionless  form, 

aT/PQ  , at  an  angle  of  inclination  of  30  degrees.  It  is  clear  from  this 

figure  that  the  radial  stress  attenuates  to  a value  of  o /P  = 0.25  at 

r o 

r = 2r  . 
o 

The  distribution  of  tangential  stress  along  the  boundary  of  the 
borehole  whose  angle  of  inclination  is  30  degrees  is  shown  in  Figure  3.3. 

Figure  3.^+  shows  the  distribution  of  radial  and  shear  stresses  along 
the  boundary  of  the  borehole  whose  angle  of  inclination  is  30  degrees. 

It  is  interesting  to  note  that  the  shear  stress,  , along  the 

boundary  of  the  borehole  is  not  zero  but  would  be  if  the  material  were 
isotropic. 


A table  of  factors  for  converting  U.  S.  customary  units  of  measurement 
to  metric  (SI)  units  is  presented  on  page  2. 
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3.6  SPATIAL  DISPLACEMENT 

DISTRIBUTION  FOR  SAMPLE  PROBLEM 

Figure  3.5  shows  a typical  result  of  the  radial  displacements  along 
the  radius  of  the  borehole  for  0=0  and  90  degrees  at  an  angle  of  in- 
clination of  30  degrees.  The  solid  line  shows  the  radial  displacement 
for  0 = 90  degrees  and  the  dashed  line  shows  the  radial  displacement 
for  0=0  degrees.  Note  that  both  the  solid  and  the  dashed  lines  would 
coincide  if  the  material  is  isotropic. 

Figure  3.6  shows  a radial  displacement  contour  in  dimensionless 

form  u /u  (r  ) , at  an  angle  of  inclination  of  30  degrees  where  u 
r r o r 

is  the  radial  displacement  at  r and  0 , and  u^(rQ)  is  the  radial 
displacement  along  the  boundary  of  the  borehole.  It  is  clear  from 
Figures  3-5  and  3.6  that  for  this  case  approximately  one  half  of  the 
borehole  volume  change  is  due  to  strains  in  the  material  within  one 
borehole-radius  o^  the  sidewall.  It  is  also  clear  that  three  fourths  of 
the  borehole  volume  change  is  due  to  strains  that  occur  at  less  than 
three  radii  from  the  borehole  sidewall.  Since  the  radius  of  a borehole 
pressuremeter  test  is  typically  1.5  inches,  only  a very  small  volume  of 
in  situ  material  close  to  the  borehole  can  significantly  influence  the 
test  results. 

The  distribution  of  the  radial  and  tangential  displacements  along 
the  boundary  of  the  borehole  whose  angle  of  inclination  is  30  degrees  is 
shown  in  Figure  3*7.  As  the  figure  shows,  the  borehole  deforms  to  an 
elliptical  shape  under  load. 

Figure  3.8  shows  a tangential  displacement  contour  in  dimensionless 
form  vQ/Vg(ro)  where  vQ  is  the  tangential  displacement  at  r and 
0 , and  Vg(rQ)  is  the  tangential  displacement  along  the  boundary  of 
the  borehole. 
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Figure  3.1  Distribution  of  radial  and  tangential  stresses 
along  the  radius  for  0=0  degrees  at  an 
angle  of  inclination  of  30  degrees . 
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Distribution  of  radial  and  shear  stresses  along  the  bourn 
borehole  whose  angle  of  inclination  is  30  degrees. 


RADIAL.  DISPLACEMENT,  INC 


RADIUS,  INCHES 

Figure  3-5  Distribution  of  radial  displacement  ur  along 
the  radius  for  0=0  and  90  degrees  at  an 
angle  of  inclination  of  30  degrees. 
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CHAPTER  1+ 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  solution  presented  herein  can  be  used  in  the  analysis  of  bore- 
hole pressuremeter  test  data  to  provide  an  appropriate  set  of  linear 
elastic  transverse-isotropic  constitutive  properties  for  a given  medium 
and  to  provide  an  index  of  a specific  site's  degree  of  anisotropy.  To 
do  this,  a series  of  pressuremeter  tests  in  boreholes  inclined  at 
several  different  angles  to  the  material's  axis  of  symmetry  have  to  be 
conducted. 

It  is  recommended  that  this  solution  be  used  at  a very  low  stress 
.evel  or  whenever  the  material  of  interest  is  assumed  to  be  linear 
elastic  transverse-isotropic.  For  a highly  nonlinear  material,  however, 
this  solution  gives  effective  constitutive  properties. 
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APPENDIX  A:  NOTATION 
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Material  properties  matrix 

Arbitrary  constants 

Radius  of  the  cylindrical  cavity 

Elastic  properties  for  the  coordinate  system  xyz 

Elastic  properties  for  the  coordinate  system  x'y'z' 

Arbitrary  constant 

Complex  operator 

Young ' s modulus  in  the  plane  of  isotropy 

Young's  modulus  in  a plane  normal  to  the  plane  of 
isotropy 

Shear  modulus  for  the  plane  of  isotropy 

Shear  modulus  for  a plane  normal  to  the  plane  of 
isotropy 

Differential  operators  of  the  second,  third,  and 
fourth  orders,  respectively 

Differential  operator  of  the  sixth  order 

Complex  stress  function 

Complex  conjugate  of  P^( z^) 

Applied  load  on  the  boundary  of  the  borehole 
Stress  function 

Transformation  matrix  for  this  appears  in 
Equation  2.3 

Complex  stress  function 

Complex  conjugate  of  (^(z^) 

Stress  function 
Cylindrical  coordinate  system 
Arc  length  along  the  boundary 
Radial  displacement 

Radial  displacement  along  the  boundary  of  the 
borehole 

Displacements  in  the  x-,  y-,  and  z-directions , 
respectively 

Initial  volume  of  a unit  length  of  the  borehole 
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Tangential  displacement 

Tangential  displacement  along  the  boundary  of  the 
borehole 

Cartesian  coordinate  system  in  which  the  solution 
of  the  borehole  pressuremeter  problem  is  solved 

Cartesian  coordinate  system  in  which  x'y’  is 
parallel  to  the  plane  of  isotropy 

The  x component  of  a distributed  surface  force 
per  unit  area 

The  y component  of  a distributed  surface  force 
per  unit  area. 

Complex  plane 
Complex  conjugate  of  z 

.K. 

Change  in  volume  per  unit  length  of  a borehole 

Total  normal  strain  components  parallel  to  x-,  y-, 
and  z-axes,  respectively 

Total  normal  strain  components  parallel  to  x'-, 
y'-,  and  z'-axes,  respectively 

Total  shearing  strain  components  in  xz- , yz-, 
and  xy-planes,  respectively 

Total  shearing  strain  components  in  x'z'-,  y'z'-, 
and  x'y' -planes,  respectively 

Root  of  the  algebraic  equation  that  corresponds  to 
differential  Equations  3.9  and  3.10 

Poisson's  ratio  that  characterizes  the  transverse 
reduction  in  the  plane  of  isotropy  due  to  stress 
in  the  same  plane 

Poisson's  ratio  that  characterizes  the  transverse 
reduction  in  the  plane  of  isotropy  due  to  stress 
normal  to  it 

Radial  stress 

Total  normal  stress  components  parallel  to  the  x-, 
y-,  and  z-axes,  respectively 

Total  normal  stress  components  parallel  to  the  x'-, 
y'-,  and  z'-axes,  respectively 

Total  shearing  stress  components  in  xz-,  yz-,  and 
xy-planes,  respectively 

Total  shearing  stress  components  in  x'z'-,  y'z'-, 
and  x'y'-planes,  respectively 
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oQ  Tangential  stress 

n Shear  stress  in  0z  plane 

0z 

iji  Angle  of  inclination  of  the  borehole 
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